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We present exact solutions of the massless Klein-Gordon equation in a spacetime in which an inh- 
nite straight cosmic string resides. The first solution represents a plane wave entering perpendicular 
to the string direction. We also present and analyze a solution with a static point-like source. In 
the short wavelength limit these solutions approach the results obtained by using the geometrical 
optics approximation: magnification occurs if the observer lies in front of the string within a strip 
of angular width SnG/j,, where p is the string tension. We hnd that when the distance from the 
observer to the string is less than 10“^(G/i)~^A ~ 150Mpc(A/AU)(G/i/10“®)“^, where A is the 
wave length, the magnification is significantly reduced compared with the estimate based on the 
geometrical optics due to the diffraction effect. For gravitational waves from neutron star(NS)-NS 
mergers the several lensing events per year may be detected by DECIGO/BBO. 

PACS numbers: 98.80.Cq, 97.60.Lf 


I. INTRODUCTION 

Typical wavelength of gravitational waves from astro- 
physical compact objects such as BH(black hole)-BH bi¬ 
naries is in some cases very long so that wave optics must 
be used instead of geometrical optics when we discuss 
gravitational lensing. More precisely, if the wavelength 
becomes comparable or longer than the Schwarzschild 
radius of the lens object, the diffraction effect becomes 
important and as a result the magnification factor ap¬ 
proaches unity U i i i i Mainly due to the 
possibility that the wave effects could be observed by 
future gravitational wave observations, several authors 

IE IE lE IE ^ 1^ IS Q Si studied wave ef¬ 

fects in gravitational lensing in recent years. 

In most of the works which studied gravitational lens¬ 
ing phenomenon in the framework of wave optics, iso¬ 
lated and normal astronomical objects such as galaxies 
are concerned as lens objects. Recently Yamamoto and 
TsunodaS studied wave effects in gravitational lensing 
by an infinite straight cosmic string. The metric around 
a cosmic string is completely different from that around 
a usual massive object. 

Cosmic strings generically arise as solitons in a 
grand unified theory and could be produced in the 
early universe as a result of symmetry breaking phase 
transition^^ ^|. If symmetry breaking occurred after 
inflation, the strings might survive until the present uni¬ 
verse. Recently, cosmic strings attract a renewed inter¬ 
est partly because a variant of their formation mecha¬ 
nism wa s propo sed in the context of the brane inflation 
scenario [TEITE1^ 1^ 1^ 1^ 1 ^ . In this scenario infla¬ 
tion is driven by the attractive force between parallel D- 
branes and parallel anti D-branes in a higher dimensional 
spacetime. When those brane-anti-brane pairs collide 
and annihilate at the end of inflation, lower-dimensional 
D-branes, which behave like monopoles, cosmic strings 
or domain walls from the view point of four-dimensional 
observers, are formed generically m EE mm EE. 


For some time, cosmic string was a candidate for 
the seed of structure formation of our universe, but 
this possibility was ruled out by the measurements of 
the spectrum of cosmic microwave background (CMB) 
anisotropies EE |^. The current upper bound on the 
dimensionless string tension Gp is around 10“^ ~ I 0~^, 
which comes from the observations of CMB EE 1^ 1^1^ 
and/or the pulsar timing EEmmEE Although cos¬ 
mic string cannot occupy dominant fraction of the energy 
density of the universe, its non -neg ligible population is 
still allowed observationallvpil l41|. In fact, Sazhin et 
al.EEEEl reported that CSL-1, which is a double image 
of elliptical galaxies with angular separation 1.9 arcsec, 
could be the first case of the gravitational lensing by a 
cosmic string with Gp « 4 x 10“^. 

We study in detail wave effects in the gravitational 
lensing by an infinite straight cosmic string. In Ref. 0: 
wave propagation around a cosmic string was studied but 
they put the waveform around the string by hand. ^ 
Their prescription is correct only in the limit of geomet¬ 
rical optics, which breaks down when the wavelength be¬ 
comes longer than a certain characteristic length. In this 
paper, we present exact solutions of the (scalar) wave 
equation in a spacetime with a cosmic string. We an¬ 
alytically show that our solutions reduce to the results 
of the geometrical optics in the short wavelength limit. 
We derive a simple analytic formula of the leading or- 


^ After submitting this paper, we have noticed a paper El in 
which the solutions of the wave equations around the cosmic 
string are given, though the apparen t ex pressions are different 
from those given in this paper. In l44l the author estimated 
the amplitude of the diffracted wave to be suppressed by 0{Gfi) 
compared with that corresponding to the geometrical optics. We 
show that the importance of the diffraction effects are determined 
by the combination of three parameters, G/r, the distance from 
the string to the observer and the wavelength and that the rel¬ 
ative amplitude of the diffracted wave can be 0(1) for realistic 
astrophysical situations. 
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der corrections to the geometrical optics due to the finite 
wavelength effects and also an expression for the long 
wavelength limit. Interference caused by the lensing re¬ 
mains due to the diffraction effects even when only a 
single image can be seen in the geometrical optics. This 
fact increases the lensing probability by cosmic strings. 

This paper is organized as follows. In section II, we 
construct a solution of the wave equation on a back¬ 
ground spacetime with an infinite straight cosmic string 
in the case that a source of the wave is located infinitely 
far. An extension to the case in which a point source is 
located at a finite distance is given in Appendix B. In sec¬ 
tion III, we study properties of the solution obtained in 
sec. II in detail. In section IV, we focus on compact bina¬ 
ries as the sources of gravitational waves and discuss the 
possible effects due to finiteness of the lifetime and the 
frequency evolution of the binaries on the detection of the 
gravitational waves which pass near a cosmic string. We 
also give a rough estimate for the event rate of the lens¬ 
ing of gravitational waves from NS-NS mergers assuming 
DECIGO/BBO. Section V is devoted to summary. 


II. A SOLUTION OF THE WAVE EQUATION 
AROUND AN INFINITE STRAIGHT COSMIC 
STRING 


A solution of Einstein equations around an infinite 
straight cosmic string to first order in G/i is given by 


d'^s = -dt^ + dr^ -b (1 - + dz^, (1) 

where (r, z, 0) is a cylindrical coordinate(0 < 9 < 27r) 
and 27rA « SttG/t is the deficit angle around the cosmic 
string. Spatial part of the above metric describes the 
Euclidean space with a wedge of angular size 27rA re¬ 
moved. Due to the deficit angle around a string, double 
images of the source are observed with an angular sepa¬ 
ration < 27rA when a source is located behind the string 
in the limit of geometrical optics. In general for a wave 
with a finite wavelength, some interference pattern ap¬ 
pears. An exact solution of Einstein equations around a 
finite thickness string has been already obtained |4^ , but 
we use the metric o as a background since the string 
thickness is negligibly small compared with the Einstein 
radius, « ttDA, where D is the distance from the ob¬ 
server to the string. 

Throughout the paper, we consider waves of a massless 
scalar field instead of gravitational waves for simplicity, 
but the wave equations are essentially the same in these 
two cases. An extension to the cosmological setup is 
straightforwardly done by adding an overall scale factor. 
In that case the time coordinate t is to be understood 
as the conformal time. The wave equation remains un¬ 
changed if we consider a conformally coupled field, but 
it is modified for the other cases due to curvature scat¬ 
tering. The correction due to curvature scattering of the 


Friedmann universe is suppressed by the square of the ra¬ 
tio between the wavelength and the Hubble length, which 
can be neglected in any situations of our interest. 

Our goal of this section is to construct a solution of the 
wave equation which corresponds to a plane wave injected 
perpendicularly to and scattered by the cosmic string. 
This situation occurs if the distance between the source 
and the string is infinitely large. In order to construct 
such a solution, we introduce a monochromatic source 
uniformly extended in the z-direction and localized in 
r — 9 plane, 

^ (2) 

where uj is the frequency and we have introduced H, a 
constant independent of A, to adjust the overall normal¬ 
ization when we later take the limit ro —> cx). The factor 
(1 — A)~^ appears because 0-coordinate used in the met¬ 
ric o differs from the usual angle 

^=(1-A)0. (3) 


Here we consider a uniformly extended source instead of 
a point source since the former is easier to handle. When 
the limit Tq —> oo is taken, the answers are identical in 
these two cases. The case with a point-like source at a 
finite distance is more complicated. This case is treated 
in Appendix B. 

Now the wave equation that we are to solve is 

V9r2 r 9r (1 - 90^ + ^ j ^ 

= Y3^<5(r-ro)(5(0-7r). (4) 

Since —9) satisfies the same equation 0 as (j){r, 9) 
does, (j)(r, 9) is even in 0. Thus, it can be expanded as 


(/i(r, 0) = ^ fm{r) cosm9. (5) 

m—O 


From Eqs. 0 and ®, the equations for fmii") are 


^ d, 2 \ r r \ 

- ^ ] frnir) 

dr‘^ r dr r 


— 6^} 


(-l)"* B 
' 1 - A ^ 


S{r - To), 


( 6 ) 


where Co = l,em = 2(m > 1) and = (1 — A) ^m. 
The solution of Eq. 0 except for r = Tq is a linear 
combination of Bessel function and Hankel function. We 
impose that the wave (p is regular at r = 0 and pure 
out-going at infinity. Further, imposing that the wave is 
continuous at r = Tq, fmir) becomes 


fm{r) = Nm {Hl^2i.'^ro)Jv^{ujr)Q{ro - r) 

- ro)j , 


( 7 ) 






where 0(a:) is the Heaviside step function. Substitut¬ 
ing Eq. o into Eq. the normalization factor is 
determined as 
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N„ 


B £^(-ir 

1 — A 2ttuj 


(coro) - 

Broe^i-1)^ 


4z(l-A) 


(wro) 

( 8 ) 


where ' denotes a differentiation with respect to the argu¬ 
ment. Erom Eqs. o and 0 with the aid of the asymp¬ 
totic formulae of the Bessel and Hankel functions, (/)(r, 6) 
for To oo can be written as 


(j){r, 9) 


-iB 



OO 

_ im.7rA 

emi^e {lot) cos mO. 

m—0 


(9) 


We determine the overall normalization of the source am¬ 
plitude B, independently of G/i, so that Eq. 0 becomes 
a plane wave when G/i = 0. This condition leads 

to i? = Then, finally (p becomes 


p{r, 9) 


1 

1 - A 


CXD 

_ ‘imA tt 

em*'"e 2(1-*) Jjy„(wr) cosm6». 

TTi — O 


( 10 ) 


III. LIMITING BEHAVIORS OF THE 
SOLUTION 

A. Approximate waveform in the wave zone 

The solution describes the waveform propagating 
around a cosmic string. But it is not easy to understand 
the behavior of the solution because it is given by a series. 
In fact, it takes much time to perform the summation in 
Eq. 1101) numerically for a realistic value of tension of the 
string, say, G^jl < 10“® because of slow convergence of the 
series. In particular it is not manifest whether the ampli¬ 
fication of the solution in the short wavelength limit co¬ 
incides with the one which is obtained by the geometrical 
optics approximation. Therefore it will be quite useful if 
one can derive a simpler analytic expression. Here we re¬ 
duce the formula by assuming that the distance between 
the string and the observer is much larger than the wave 
length, 


^ = wr > 1, 


( 11 ) 


which is valid in almost all interesting cases. 

Using an integral representation of the Bessel function. 




(12) 


where the contour of the integral G is such as shown in 
Eig.0 Eq. 10 can be written as 


CXD 

X ^ -f e-*’"®).(13) 

m—0 


When t is in the segment of the integration contour G 
along the imaginary axis, the summation over m does 
not converge because the absolute value of each term in 
the summation is all unity. In order to make the se¬ 
ries to converge, we need to think that the integration 
contour G is not exactly on the imaginary axis but t al¬ 
ways has a positive real part. Eor bookkeeping purpose, 
we multiply each term in the sum by a factor (e 

is an infinitesimally small positive real number). Then 
Eq. itT^ becomes 


•Pit 0) = 0) + V^(^, -9), 


where 9) is defined by 


me) 


1 1 /■ , e^sinht 


with 





a{9) 


(14) 


(15) 


(16) 


and a{9) := (ttA — (1 — A)9)^/^. 

Now we find that all we need to evaluate is V'(C)^) I'd 
order to obtain an approximate formula for 9). This 
integral will not be expressed by simple known functions 
in general, but the integration can be performed by using 
the method of steepest descent in the limit ^ ^ I. 

The integrand of Eq. OSl) has two saddle points located 
at t = t_|_ z7r/2 and t = « —i7r/2 in the vicinity of 

the integration contour G. We should also notice that 
the integrand has a pole at t = t*, which is also infinites¬ 
imally close to the contour of the integral G. This pole 
is located near the saddle point at t as far as A 

and 9 are small. Hence the treatment of the saddle point 
at t = is much more delicate than that of the saddle 
point at t = t_. We only discuss the saddle point at 
t = then the case at t = is a trivial extension. 

When 3fi(t) > 0,3(t) < or ^(t) < 0,3(t) > i-|, 
which corresponds to shaded regions in Fig. □ g€sinht 
verges in the limit ^ ^ oo. If a{9) > 0, the pole at t = t* 
is in the bottom-left unshaded region. In this case we 
cannot deform the contour to the direction of the steep¬ 
est descent at t = without crossing the pole at t = G. 
The deformed contour which is convenient to apply the 
method of the steepest descent is such that is shown as 
G in Fig. ^ When we deform the integration contour 
from G to G, there arises an additional contribution cor¬ 
responding to the residue at t = t* when a{9) > 0. On 
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the other hand, if a(0) < 0, the pole is in the top-left 
shaded region. In this case, we can deform the contour 
of the integral to the direction of the steepest descent 
without crossing the pole t*. Hence no additional term 
arises. 

From these observations, we find that it is necessary 
to evaluate the integral 03 separately depending on the 
signature of a{d). Though the calculation itself can be 
done straightforwardly, it is somewhat complicated be¬ 
cause the saddle point and the pole are close to each 
other. When the pole is located inside the region around 
the saddle point that contributes dominantly to the in¬ 
tegral, a simple Gaussian integral does not give a good 
approximation. Detailed discussions about this point are 
given in Appendix A. Here we only quote the final result 
which keeps terms up to 0(l/v^), 




where 

<5(0) 

ct ( 0 ) 

and 


exp cos ®(q^(^)) 

cr(0) / ia(d) i 


xErfc( 


-b 


1 1 




V2< 1 - A 1 _ gT^G-a(G/Vf) ^ 


(17) 


*(1-A)v^ 

sign(a(0)). 


1 — exp ( i 


a(0) 


(l-A)V?;j 


,( 18 ) 

(19) 


r*+oo 


Erfc(a;) ■= dte 


( 20 ) 


We are mostly interested in the cases with A, 9 1. 

Then, we have a{9)/^ ^ 1, and therefore <5(0) reduces 
to a(0). The second term in Eq. If 711 is the contribution 
from the integral around the saddle point at t = t_|_ along 
the contour C. This term is not manifestly suppressed by 
1/v^. As far as a{9) is fixed, this term does not vanish 
in the limit ^ ^ oo. Of course, if we fix A and 0 first, and 
take the limit ^ oo, the argument of the error function 
goes to -boo and the function itself vanishes. However, 
a(0) vanishes at 0 = 7rA/(l — A). Hence even for a very 
large value of ^ there is always a region of 0 in which 
this second term cannot be neglected. However, for 0 
in such a region, a(0) cannot be very large. Therefore, 
we can safely drop the second term in the exponent. On 
the other hand, the last term in Eq. 03, which is the 
contribution from the saddle point at t = t_, is always 
suppressed by l/\/^. Hence, this term does not give any 
significant contribution for ^ ^ 1. The first term in 
Eq. (I14II can be dropped in the same manner for ^ ^ 1. 
Keeping only the terms which possibly remain in the limit 
^ > oo, we finally obtain 




exp 



0(a(0)) 


^g*«-5«"(0Erfc 

VTT 

+ ( 0 -- 0 ). 


H^-in/4 

V2 

( 21 ) 


For illustrative purpose, we compared the estimate 
given in Eq. m with the exact solution Eq. m in 
Fig. El They agree quite well at ^ » 1. The deficit angle 
and the observer’s direction are chosen to be A = 0.0025 
and 0 = 0, respectively. 



FIG. 1: Black, dotted and dashed lines are contours of the 
integral C, C and Ch, respectively. ±i^ are the saddle points 


of 



FIG. 2: Configuration of the source, the cosmic string and 
the observer. A and B are the positions of a source. O and 
P are the positions of the cosmic string and the observer, 
respectively. In this figure, the wedge AOB is removed and 
thus A and B must be identified. 


B. Geometrical optics limit 

Geometrical optics limit corresponds to the limit ^ ^ 
oo with A and 0 fixed. In this limit a{9) also goes to -boo. 
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€ 0 

FIG. 3: Comparison between the exact solution Eq. GOJ and the approximate one Eq. I17II . ttA is 0.0025. Black line and 
dotted one correspond to the exact solution and the approximate one, respectively. We see that except for small ^ the dotted 
line overlaps the black one. In the right panel, the relative error is about 10“®. 


and hence the error function in Eq. cu vanishes. Hence 
the waveform in the geometrical optics limit, which we 
denote as (^go, becomes 

( 22 ) 

where ip is defined by Eq. 0 . 

Since 4> and hence ^go are even in 9, it is sufficient to 
consider the case with 0 > 0. In Eig.|21 the configuration 
of the source, the lens and the observer is drawn in the 
coordinates in which the deficit angle 27rA is manifest, 
i.e., the wedge ADB is removed from the spacetime. Both 
points A and B indicate the location of the source. The 
fines DA and DB are to be identified. The angle made 
by these two lines is the deficit angle. The locations of 
the string and the observer are represented by 0 and P, 
respectively. In our current setup the distance between 
0 and A (= Tq) is taken to be infinite. When p > ttA, 
only the source A can be seen from the observer. This 
corresponds to the fact that only the first term remains 
for p > ttA in Eq. 123). Eor p > ttA, we have 

(23) 

This is a plane wave whose traveling direction is p = 
—ttA, which is the direction of AP in Eig. |21in the limit 
To = |A0| ^ oo. 

For \p\ < ttA, (^go is 

{i,0) = . (24) 

This is the superposition of two plane waves whose trav¬ 
eling directions are different by the deficit angle 27rA. 
Hence amplification of the images and interference occur 
for \p\ < ttA as expected. 

As we shall explain below, Eq. (|23 coincides with the 
one derived under the geometrical optics. In geometrical 
optics, wave form is given by El 

<pgo = eyi'p[iu}T{xj) — i'Krij], (25) 

j 


where x represents a two-dimensional vector on the lens 
plane and T{x) represents the summation of time of flight 
of the light ray from the source to the point x on the lens 
plane and that from the point x to the observer. Xj 
is a stationary point of T{x), and nj = 0,1/2,1 when 
Xj is a minimum, saddle and maximum point of T{x), 

respectively. The amplitude ratio |u(x)| ' is written as 

u{x) = 1/ det[5ab - dadbipix)], (26) 

where ipix) in Eq. (I26II is the deflection potential ^3 
which is the integral of the gravitational potential of the 
lens along the trajectory between the source and the ob¬ 
server. Eq. 1251) represents that the wave form is obtained 
by taking the sum of the amplitude ratio |m(x,)| ' of 
each images with the phase factor ^ 

lens is the straight string, the spacetime is locally flat 
everywhere except for right on the string. This means 
that the deflection potential 4’i^) is zero and hence the 
amplitude ratio is unity for all images and the tra¬ 
jectory where the time of flight T{x) takes the extremal 
value is a geodesic in the conical space, and T(x) of any 
geodesic takes minimum, which means rij = 0. There are 
two geodesics if the observer is in the shaded region in 
Fig. El The time of flight along the trajectory AP is 

Ta = lim 1^1 ro-I-r cos(7rA-I-(/j), (27) 

ro—*oo 

where r = |0P|. The time of flight along the trajectory 
BP is obtained by just replacing p with —p. Hence, sub¬ 
stituting m into E3), we find that the waveform in 
the geometrical optics is the same as Eq. EH) except for 
an overall phase This factor has been already ab¬ 

sorbed in the choice of the normalization factor B in our 
formula m- 

We define the amplification factor 

4>VL {^,0)’ 




( 28 ) 
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where ipvh is the unlensed waveform. Using Eq. it^ . the 
amplification factor of (j)go for |</3| < ttA is given by 

Fgo(?, 0) « cos(7rAe(^), (29) 

where we have assumed (p and A are small and dropped 
terms higher than quadratic order. It might be more 
suggestive to rewrite the above formula into 

l^go(C.6')l « 2cos(7rAw2/), (30) 

where y = r sin p. The distance from a node to the next 
of when the observer is moved in y-direction is A/ttA, 
where A is a wavelength. This oscillation is seen in the 
right panel of Fig. 0 


C. Quasi-geometrical optics approximation 

In the previous subsection, we have derived the wave¬ 
form in the limit |a(±d)| —> oo which corresponds to 
the geometrical optics approximation. Here we expand 
the waveform laj to the lowest order in l/a(±0). This 
includes the leading order corrections to the geometrical 
optics approximation due to the finite wavelength effects. 

For the same reason as we explained in the previous 
subsection, we assume that A and p are small. Using 
the asymptotic formula for the error function Eq. 
the leading order correction due to the finite wavelength, 
which we denote as (5^!)qgo, is obtained as 

Hq,go{i,0) = 


p^+iir/A 


1 


-b 


\oi{0) 

e*^+”/4 27rA 

(7rA)2 — (^2 ’ 


,[- 0 ) 


(31) 


T{x), which we call the image points, contribute to the 
waveform. The first length scale is = |a(±0)|/v^ x r 
which is defined as the separation between an image point 
and the string on the lens plane. In this picture we expect 
that paths whose pathlength is longer or shorter than the 
value at an image point by about one wavelength will not 
give a significant contribution because of the phase can¬ 
cellation. Namely, only the paths which pass within a 
certain radius from an image point need to be taken into 
account. Then such a radius will be given by rp = -s/Ar, 
which we call Fresnel radius. Namely, a wave with a fi¬ 
nite wavelength can be recognized as an extended beam 
whose transverse size is given by rp. The ratio of these 
two scales gives a{±9)-. 


|a(±6l)| 


V^rs 

rp 


When Xg ^ rp , i.e., a{±6) » 1, the beam width is 
smaller than the separation. In this case the beam image 
is not shadowed by the string, and therefore the geo¬ 
metrical optics becomes a good approximation. When 
Ts ^ rp, i.e., 


a{±0) < 1, (33) 

we cannot see the whole image of the beam, truncated 
at the location of the string. Then the diffraction ef¬ 
fect becomes important. The ratio of the beam image 
eclipsed by the string determines the phase shift and the 
amplification of the wave coming from each image. If we 
substitute |(/3| ~ 0 as a typical value, we obtain a rough 
criterion that the diffraction effect becomes important 
when 


A ^ 27r(7rA)^r, 


(34) 


As is expected, the correction blows up for \p\ « ttA, 
where a(0) or a{—6) vanishes, irrespectively of the value 
of In such cases, we have to evaluate the error function 
directly, going back to Eq. m- 

The expression on the first line in Eq. m manifestly 
depends only on a{zt9) aside from the common phase 
factor e*^. This feature remains true even if we consider 
a small value of a(±0). This can be seen by rewriting 
Eq. ll^ as 

(®) /— a(0)\ 

<P{^,0) ^ ^ Erfc +{0^ -0)m 

The common phase does not affect the absolute mag¬ 
nitude of the wave. Except for this unimportant overall 
phase, the waveform is completely determined by a{±0). 

The geometrical meaning of these parameters a{±9) is 
the ratio of two length scales defined on the lens plane. 
To explain this, let us take the picture that a wave is 
composed of a superposition of waves which go through 
various points on the lens plane. In the geometrical op¬ 
tics limit the paths passing through stationary points of 


or ■C ^ (^^)~^ fo terms of 

The same logic applies for a usual compact lens object. 
In this case the Fresnel radius does not change but the 
typical separation of the image from the lens is given 
by the Einstein radius rp ^ \/AGMr, where M is the 
mass of the lens. Then the ratio between rp and rp is 
given by rpjrp = ^jGMjX, which leads to the usual 
criterion that the diffraction effect becomes important 
when A>GM 0,1113,1113. 

From the above formula m, we can read that the 
leading order corrections scales like oc \/Xlr. This de¬ 
pendencies on A and r differ from the cases that the lens 
is composed of a normal localized object, in which the 
leading order correction due to the finite wavelength is 
0{X/M) [13. 

The condition for the diffraction effect to be impor¬ 
tant II33II can be also derived directly from Eq. 1131 |l . In 
order that the current expansion is a good approxima¬ 
tion, (j^qgo must be smaller than (jjgo- This requires that 
l/a(±d) 1, which is identical to ijSSI)- 

We plot the absolute value of the amplification fac¬ 
tor under the quasi-geometrical optics approximation as 
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FIG. 4: The absolute value of the amplification factor as a 
function of ^ for 6 — 0. Black line and dashed one correspond 
to Eq. 12111 and the quasi-geometrical optics approximation, 
respectively. The string tension is chosen to be Gfi = 10~^. 


dashed line in Fig.0| We find that the quasi-geometrical 
optics approximation is a good approximation for ^ ^ 
A“^. For ^ the quasi-geometrical optics approx¬ 

imation gives a larger amplification factor than the exact 
one. 

In the quasi-geometrical optics approximation, we find 
from Eqs. ll^ and the absolute value of the ampli¬ 
fication factor for </? = 0 is 


| F (^, 0)|«2 

From this expression, we find that the position of the first 
peak of the amplification factor lies at ^ ~ 4.25 x (ttA)”^, 
which can be also verified from Fig. IS For ^ < A ^ the 
present approximation is not valid, but we know that the 
amplification factor should converge to unity in the limit 
^ > 0, where rp is much larger than rg. 

We show in Fig. 0 the absolute value of the amplifica¬ 
tion factor as a function of ip for four cases of ^ around 
A“^. Top left, top right, bottom left and bottom right 
panels correspond to ^(ttA)^ = 0.5,1,2 and 4, respec¬ 
tively. Black curves are plots for Eq. m and the dotted 
ones are plots for the quasi-geometrical optics approxi¬ 
mation. As is expected, the error of the quasi-geometrical 
optics approximation becomes very large near ip = ttA, 
where a{6) vanishes. As the value of ^ increases, the an¬ 
gular region in which the quasi-geometrical optics breaks 
down is reduced. 

Interestingly, the absolute value of the amplification 
factor deviates from unity even for ip > ttA which is 
not observed in the geometrical optics limit. This is a 
consequence of diffraction of waves, the amplitude of os¬ 
cillation of the interference pattern becomes smaller as 
9 becomes larger, which is a typical diffraction pattern 
formed when a wave passes through a single slit. The 
broadening of the interference pattern due to the diffrac¬ 
tion effect means that the observers even in the region 
|(p| > ttA can detect signatures of the presence of a cos¬ 
mic string. 


L - \ ---2 + T 

y TT^inA^ V2 4 


1/2 




But the deviation of the amplification from unity out¬ 
side the wedge ip > ttA is rather small except for the 
special case ^(ttA)^ I: for ^(ttA)^ <C 1 the magnifica¬ 
tion is inefficient and for ^ 1 the magnification 

itself does not occur. Hence the increase of the event 
rates of lensing by cosmic strings compared with the es¬ 
timate under the geometrical optics approximation could 
be important only when the relation ^(ttA) ft! 1 is sat¬ 
isfied. If we take D = 10^®cm and uj = 10“^Hz which 
is in the frequency band of LISA(Laser Interferometer 
Space Antenna) , we find that the typical value of Gp, 
is ft 2 X 10-9. 

So far, we have considered the stringy source rather 
than a point source. Extension to a point source can be 
done in a similar manner to the case of the stringy source 
and is treated in Appendix B. The result is 


(36) 

where D = \/{f + Tq)^ + is the distance between 
the source and the observer. T, which is defined by 
Eq. fgUll . is related to ip as 


T{x,9)=[e-^^i;{^,9)] + {9 ^-9). 


(37) 


Hence p for the point source is similar to that for the 
stringy source. In particular, assuming that A,(/? <C 1, 
and keeping terms which could remain for tor, coro S> 1, 
we have 


D. Simpler derivation of Eq. ll^ . 

We have derived an approximate waveform which 
is valid in the wave zone from the exact solution of the 
wave equation Eq. m- Here we show that Eq. can 
be obtained by a more intuitive and simpler method. In 
the path integral formalism 0. the wave form is given 
by the sum of the amplitude exp {iujT^s)) for all possible 
paths which connect the source and the observer. Here 
T(s) is the time of flight along the path s. If the cosmic 
string resides between the source and the observer, the 
wave form will be given by the sum of two terms one 
of which is obtained by the path integral over the paths 
which pass through the upper side of the string {y > 0) in 
Fig. 12 and the other through the lower side of it (y < 0). 
The waveform coming from the former contribution will 
be given by 


/ OO pOO _^ 

dzQ / d2/qe*“6AQ|+|qP|)^ (39) 

-00 J 0 
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FIG. 5: Black line and dotted one correspond to Eq. (HU and the quasi-geometrical optics approximation, respectively. The 
string tension is chosen to be Gfi — 10~®. 


where Q = (0, Uq, zq) is a point on the lens plane specified 
by X = 0. One can determine the normalization constant 
A by a little more detailed analysis, but we do not pursue 
it further here. By integrating Eq. we recover the 
first term in Eq. (ESl. 


E. Long wavelength limit 

For completeness, we consider the case in which the 
wavelength is longer than the distance from the string 
^ ^ 1. In this limit, the first few terms in Eq. m 
dominate, and we find 

1 A log 2—7rA/22 

~ — r(l +A) —^^^^cos6l. (40) 

In particular, for ^ > 0 Eq. (Eni becomes (1 — A) ^ 

which is larger than unity. This differs from the cases of 
gravitational lensing by a normal compact object, where 
the amplification becomes unity in the long wavelength 
limit. The reason why the amplification differs from unity 
even in the long wavelength limit is that the space has 
a deficit angle and hence the structure at the spatial in¬ 
finity is different from the usual Euclidean space. Waves 


with very long wavelengths do not feel the local struc¬ 
ture of string. However, uniform amplification of waves 
should occur as a result of total energy flux conservation 
because the area of the asymptotic region at a constant 
distance from the source is reduced due to the deficit 
angle. In this sense such modes feel the existence of a 
string. 


IV. CONNECTIONS TO OBSERVATIONS 
A. Compact binary as a source 

In this section, we consider compact binaries as sources 
of gravitational waves. Gravitational waves from com¬ 
pact binaries are clean in the sense that the waves are 
almost monochromatic: the time scale for the frequency 
to change is much longer than the orbital period of the 
binary except for the phase just before plunge. Hence in¬ 
terference between two waves coming from both sides of 
the cosmic string could be observed by future detectors. 

Since each compact binary has a finite lifetime, lensing 
events can be classified roughly into two cases. If the dif¬ 
ference between the times of flight along two geodesics is 
larger than the lifetime of the binary, we will observe two 
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independent waves separately at different times. On the 
other hand, if the time delay is shorter than the lifetime, 
what we observe is the superposition of two waves. 

The remaining lifetime of the binary Tnfe when the pe¬ 
riod of the gravitational waves measured by an observer 
is Pgw is estimated as 


where we assumed that S{t) is real. Substituting Eq. ll^ 
to the above expression, we have 

+C.C. (46) 


Tlife 


9.2 X 10 


-4 


(1 + v) 


1/3 


(1 + V 



5/3 

Pgw, 


. 

where 77 is the mass ratio of the binary (ry < 1 ), M is the 
mass of the more massive star in the binary and zs is the 
source redshift. 

The time delay Tdeiay is 


TV 

Tdeiay « 2-^(7J7rA. (42) 


Taking the typical values of parameters as rVo/D — IGpc 
and ip = ttA, the condition Tnfe Tdeiay gives the upper 
bound on the mass M, 


M < 8 X 10^ 


yy3/5 



Pgw \ 
lO^sec ) 


8/5 

Mq. 


The time scale for the orbital frequency of the binary 
to change is the same order as Tnfe. Hence the condi¬ 
tion Tiife Tdeiay implies that the frequencies of two 
waves are almost the same. The left and right panels 
in Fig. 1^ which correspond to different frequencies of 
gravitational waves show the region where the condition 
Eq. MHII is satisfied for three different values of string 
parameter A. The shaded area represent the parameter 
region beyond the detector’s sensitivities. In the left and 
right panels we assumed, respectively, that the thresh¬ 
old value for detection in strain amplitude for LISA and 
DECIGO(DECihertz Interferometer Gravitational wave 
Observatory) ^/BBO (Big Bang Observer l[50|. which 
are given by and We find 


that both cases Tufe T/jeiay and Tufe T/jeiay can oc¬ 
cur both for LISA and BBO/DECIGO. 


B. Waveform 


Eq. 14611 is a general formula which applies to any time 
dependent source. Here we consider the special case in 
which S{t) takes the form 

S(t) = cos dt' 0(t')^ , (47) 


with 


0(t) =u)oP uJot, (48) 

where <C I and Wq > 0 are assumed. This rep¬ 

resents a quasi-monochromatic source with its frequency 
slowly changing. Then is 


Su, = 


y/2TTUJo 


.-(id+idiX 


+27r/4 




— iTrjA 


■ (49) 


Substituting Eqs. (EHll and ll^ into Eq. ( 1 ^ . and using 
the method of the steepest descent, we have 


x) 


1 _Erfc ( - ttA _ / n{T{ip))rro 


with 


+(</> ^ -/?), 


T{p)=t-D+^^{nA-pf. 


D 


(50) 

(51) 


This represents a superposition of two waves coming from 
both sides of the string whose arrival times differ by 
\T{(p) — T{—'p)\ = 2 ^ 7 rA|(/ 3 |. In the preceding subsec¬ 
tions, we study the waveforms observed in two cases with 
Tiife ^ Tdeiay and Tiife Tdeiay- 


We can easily extend our waveform m to the case 
that the frequency of the source changes in time. Let us 
write the source as j^S{t)S{r — ro)6{9 — Tr)6{z). The 
Fourier transformation of S{t) is defined by 



Denoting the solution (/)(f, x) for a monochromatic source 
obtained in the previous sections by (x), (p can be writ¬ 
ten as 

poo 

(j){t,x) = duj +C.C., (45) 

^0 


1 . Tiife ^ Tdeiay 

As we have explained in the preceding subsection, 
what we observe is a superposition of two waves in this 
case. Because the relative phase difference of these waves 
slowly increases or decreases in time due to the frequency 
change of the binary source and the optical path differ¬ 
ence between two geodesics, we will observe the beat if 
the amplitude of the integrated relative phase difference 
over observation time is larger than 0 ( 1 ). 

The condition that the beat is observed can be derived 
as follows. If we denote the total observation period by 
Tobs, then from Eq. the integrated relative phase 
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FIG. 6: Plots of regions where Eq. (El is satisfied for three different values of the string parameter. Left and right panels are 
for 10“®Hz and O.lHz which are the frequency bands LISA and DECIGO have best sensitivities. Shaded regions are plotted 
under the assumptions that the signals satisfy SN > 10, the redshift of the source is 1 and 3 year observations. 


difference is 27:l^ipDCjoTohs^ where both r and are as¬ 
sumed to be 0{D). Hence we can observe the beat if 


Tohs ^ 


1 

27: l^ifDCjo 


(52) 


Because Tufe is roughly the same as the time scale for 
the frequency of the binary to change, i.e. T^fe ~ Wo/d>o, 
Eq. lEU can be written as 


'Eobs ^ 


?life 

27: AipDuJo 


(53) 


If Tobs is fixed, e.g. Tobs 3yr for LISA, Eq. (ICT is 
written as an lower bound on M. For Fobs = 3yr and 
Pgw = lO^sec, Eq. ( 13 ^ becomes 


M > 2.6 X 


+ f t:A 

7^3/5 


Pgw \ 
lO^sec J 


13/5 

Mq. 


(54) 

We show in Fig. [7|the region where Eq. isi is satished for 
LISA with Tobs = 3yr. We find that if G/r < 2.8 x 10“® 
which is about one order of magnitude below the current 
upper bound, LISA will detect the beat of gravitational 
waves for all observable ranges in (/i, M) space as long as 


i^ife ^ ifdely 


2 . Tlife <C Tdelay 

If Tdelay ^ Tlife, we observe the waveform of either the 
first term or the second one in Eq. at a given time. 


^ Since the lensing probability is not expected to be high, we need 
a large number of events to detect a lensing event. In such a 
situation, what gravitational wave detectors can detect is a su¬ 
perposition of various waves. Hence, signal will almost always 
have beat even if we ignore the lensing effect. 



FIG. 7: Plot of the region where Eq. il3H is satisfied. The 
frequency of the gravitational waves is assumed to be 10“®Hz. 


We show in Hg. |H1 the amplihcation of the wave corre¬ 
sponding to the first term in Eq. it3ni i as a function of 
ip — Att normalized by 1 / yAjrro/D, which is nothing but 
—a{9) in the case discussed in Sec lIIIl We find that the 
amplihcation approaches zero more slowly for i/i— ttA > 0 
and oscillates around unity for ip — ttA < 0 and the angu¬ 
lar size in which non-trivial oscillations du e to the d iffrac- 
tion effect can be observed is given by 1/ \furn\jD. Since 
Tdelay ~ {tTo/D) p7:A < Tnfe <C in the present case, 
we have (rvo/D){7:A)^ > (rro/D)7:Ap ^ 1. Therefore 
this angular size of oscillation is much smaller than ttA. 
Hence it will be very difficult to detect a lensing event in 
which this diffraction effect is relevant. 


C. Estimation of the event rate 

In this section, we estimate the detection rate of the 
gravitational lensing caused by cosmic strings for planned 
gravitational wave detectors such as LISA, DECIGO and 
BBO. 

It is well known that string network obeys the scal¬ 
ing solution where the appearance of the string network 
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FIG. 8: The absolute value of the amplification factor for Tufe 
to ^ = 0.2(7rA)~^ and (ttA)”^, respectively. 

at any time looks alike if it is scaled by the horizon size. 
There are a few dozen strings spread crossing the horizon 
volume and a number of string loops . Since 

the horizon scale increases in the comoving coordinates 
as time goes, the number of strings increase if there is 
no interaction between them. However, since strings are 
typically moving at a relativistic speed, they frequently 
intersect with each other. As a result reconnection be¬ 
tween strings occurs, reducing the number of long strings 
which extend over the horizon scale. During the process 
of reduction of the number of long strings a large number 
of string loops are formed, but they shrink and decay via 
gravitational radiation. Due to the balance of two effects, 
the number of long strings in a horizon volume remains 
almost constant in time. 

The reconnection probability p is essentially 1 for 
gauge theory solitons because reconnection allows 
the flux inside the string to take an energetically favor¬ 
able shortcut. For F-strings, the reconnection is a quan¬ 
tum process and its probability is roughly estimated as 
p where Qb is the string coupling and is predicted 

in |5fi| that 

10”^<P<1. (55) 

For D-strings, the reconnection probability might be 
0-1 ^ P ^ 1 If the reconnection probability is less 
than 1, the number of long strings is expected to be p~^ 
times larger than that in the case with p = 1. Therefore 
it is expected that in the context of cosmic strings mo¬ 
tivated by superstring theory the number of long strings 
in a horizon volume can be 10^ or more. 

To estimate the event rate for the gravitational leas¬ 
ing, here we consider a compact binary (such as binary 
neutron stars and/or black holes) as a source of gravi¬ 
tational waves. There are large uncertainties about the 
event rate of MBH (massive black hole) merger detected 
by LISA or DECIGO/BBO. Several authors 
employed a model in which MBH mergers are associated 
with the mergers of host dark matter halos to estimate 
the event rate of MBH-MBH mergers. In this model, the 


Tdeiay as a function of ip. Left and right panels correspond 


event rate is dominated by halos with the minimum mass 
Mmin above which halos have a central MBH and some 
scenario predict that the event rate could reach ^ 10"^ 
events/yr. For DECIGO/BBO, the binary neuron stars 
will be observed ~ 10® events/yr. 

The probability of lensing for a single source by an infi¬ 
nite straight cosmic string both at cosmological distances 
is 

= 3 X If” ( 5 ^) ■ (56) 

Eq. (1^ is derived under the geometrical optics approxi¬ 
mation. In section HI, we found that the signal of lensing 
by cosmic strings (the interference pattern of gravita¬ 
tional waves at detectors) extends over an angular scales 
larger than the deficit angle 27rA when the diffraction ef¬ 
fect is marginally important. This is a well known fact for 
the gravitational lensing by usual stellar objects 
As we estimated in section HI, the critical distance 
below which the diffraction effect becomes important is 

™ {^) (sr) 

Therefore the probability of lensing by cosmic strings 
may be enhanced due to the diffraction effect for ttA ft! 
10“”^ at LISA band {Paw ~ 10®sec) and for ttA ft I0“® 
at DECIGO/BBO band {Paw « lOsec). 

Assuming the prospective values of the parameters 
that determines the rate of lensing events h, we obtain 



where /(> 1) denotes the numerical factor arising from 
the enhancement of the lensing probability due to the 
diffraction effect, hs = 10® is almost upper bound on 
the total event rate of neutron star mergers detectable 
by DECIGO/BBO. If the event rate is even higher, the 
number of events becomes comparable to or larger than 









12 


the number of frequency bins. Then we will not be able 
to distinguish each event, and undistinguishable signals 
become confusion noise. In the case of LISA, this bound 
on fis is even lower. Unfortunately, a large number of 
lensing events by cosmic strings can be expected only for 
marginally large 7rA(ft! G/i/4) with a small reconnection 
probability p. 

Finally we briefly comment on the validity of the as¬ 
sumption that most of cosmic strings can be treated as 
straight ones in studying gravitational lensing by them. 
In geometrical optics approximation, only light paths 
which satisfy the Fermat’s principle contribute to the am¬ 
plification factor. If we take into account the finiteness 
of the wavelength, the trajectories whose optical path 
differences are less than a few times of its wavelength 
will dominantly contribute to the amplification factor. 
In terms of the distance on the lens plane {x = 0-plane 
in Fig. 121), the optical paths within ^ V \D from the in¬ 
tersection of the geodesic will give a dominant part of the 
amplification factor. 

In the standard literature, the typical size of small- 
scale structure of a long string is given by the gravita¬ 
tional back-reaction scale ~ 50G/rt, where t is a cosmic 
time j^. But this is not an established argument and 
some recent studies suggest that the smallest size of the 
wiggles could be much smaller than 50G^t [o^, Ell • If we 
assume here that the smallest size of the wiggles is 50G/it, 
then the condition that the straight string approximation 
is good is V \D < 50G^t. Substituting the appropriate 
values of the parameters, it gives the condition, 


1 > 


bOGfxt 


= 8 X 10"^ 


7rA\"Y ^ ^ Vr 1 

Yoi3cm; 1026cm ) 


Hence approximating a cosmic string by a straight one is 
good for wide range of possible values of the parameters. 


other length scale rp, which is called Fresnel radius, is 
the geometrical mean of the wavelength and the typical 
separation among the source, the lens and the observer. 
The waveform is characterized by the ratios between 
and rp. li rp > rg, the diffraction effect becomes im¬ 
portant and the interference patterns are formed. Even 
when the image in the geometrical optics is not directly 
seen by the observer, the interference patterns remain. 
In contrast, in the geometrical optics magnification and 
interference occur only when the observer can see two im¬ 
ages which travel both sides of the string. Namely, the 
angular range where lensing signals exist is broadened by 
the diffraction effect. This broadening may increase the 
lensing probability by an order of magnitude compared 
with that estimated by using the geometrical optics when 
the distance to the source is around the critical distance 
Dc given in Eq. m- 

We finally estimated the rate of lensing events which 
can be detected by LISA and DECIGO/BBO assuming 
BH-BH or NS-NS mergers as a source of gravitational 
waves. For possible values of the parameters that de¬ 
termines the event rate such as string reconnection rate, 
string tension and the event rate of the unlensed mergers, 
the lensing event rate could reach several per yr. 
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APPENDIX A: DERIVATION OF EQ. (HTJ) 


V. SUMMARY 

We have constructed a solution of the Klein-Gordon 
equation for a massless scalar field in the flat spacetime 
with a deficit angle 27rA « SttG/t caused by an infinite 
straight cosmic string. We showed analytically that the 
solution in the short wavelength limit reduces to the geo¬ 
metrical optics limit. We have also derived the correction 
to the amplification factor obtained in the geometrical 
optics approximation due to the finite wavelength effect 
and the expression in the long wavelength limit. 

The waveform is characterized by a ratio of two differ¬ 
ent length scales. One length scale is defined as the 
separation between the image position on the lens plane 
in the geometrical optics and the string. We have two 
since there are two images corresponding to which side 
of the string the ray travels. (When the image cannot be 
seen directly, we assign a negative number to rg.) The 


Here we derive a formula Eq. from the integral 
representation of the solution Eqs. m and (Cni). As we 
explained in the sec.HI, we have to calculate the integral 
for a{9) >0 and a(0) < 0 separately. 

a. a{d) < 0 

In this case, there is no contributions from the pole t*. 
Since ^ S> 1, the integral 

Ilf * 

i_e-T^h-*fUo(G/V?+.)’ 

(Al) 

is dominated from the two regions |t±if | ^ I/-\/C) where 
is the saddle points of 

Let us first calculate the integral around iTr/2. We 
cannot apply the method of steepest descent where the 
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denominator of the integrated function is replaced with 
the value at t = because the pole t* of the inte¬ 
grand can lie in the region |t* — i7r/2| < 1/VJ and the 
denominator is no longer constant around i7r/2. 

Fortunately the integral can be approximated written 
by the special function which can be evaluated easily. 
We first do the transformation of variable such that t — 
iix jl = (u:real number) which corresponds to the 

deformation of the contour of the integral from C to C 
as shown in Fig. ^ Expanding around i7r/2 to 

second order in u and the denominator of the integral to 
the first order in u gives the integral 

1 ia{e) \ f°° ^ e-T 

s; + (T3ij^j + 

(A2) 

where a{9) is defined by Eq. itT^ . Hence we need to 
evaluate the integral 


as a constant and apply the usual saddle point method. 
This gives 


I I 

1 - a 1 _ gTA5:(A--a(6»)/V?) ■ 


(AS) 


The sum of Eqs. m and lASII gives 9) for a{9) < 0. 


b. a{d) > 0 

In this case, there is a contribution from the pole t*. 
Hence the integral is divided into the integral around the 
pole and the one whose circuit of integration is C. 

The integral around the pole gives 

(A9) 


exp cos 




I{x) 



(A3) 


where e was introduced to remember that the imaginary 
part of a; = —e*'^/'*d(0) is positive when a{9) < 0. This 
integral is given by an error function as 


I{x) = 2iy/TTe ^ /^Erfc ^ ' 

(A4) 

This can be derived by solving a differential equation 

—I(x)=—V^—xI(x), (As) 

ax 

which follows from the definition of I(x) , with the bound¬ 
ary condition that 1(0) = in. 

Asymptotic formulas for the error function are 


Erfc( 2 ;) = e 


1 1 
2z 4.z^ 

TT ^ 1 I I 

(for ~ ^ < argz < — and p| ^ oo), 


= yJTT + e 


1 1 

^ ~ 4^ 


(for — < argz < — and l^j ^ oo), 

•v/tt 1 o 

+ ^-z+-z^ + --- , 

(for |z| < 1). 


(A6) 


Using Eq. K3ll . the integral Eq. becomes 


1 

^exp 



ia{9) 

(1-A)V? 



Erfc 



(A7) 

Next let us calculate the integral Eq. around 

—i7r/2. Since the pole t* is far from —inj2, we can 
approximate the denominator of the integrated function 


The integral around z7r/2 along the trajectory C is the 
same as for a{9) < 0 and is given by Eq. (HU. The 
only difference is the signature of d(9). By changing the 
integration variable from u to —u, a:(9) is replaced with 
—d{9) and the overall signature flips. As a result we find 
that the integration along C gives 


xErfc 


iC + 




(1-A)Ve 


V2 


,3i7r/4 


)■ 



Integral Eq. EH around —in 12 is also given by Eq. lASII . 

Combining the results of subsections a. and 6., adding 
the similar terms and also using the asymptotic 

form of Jo(^), we have Eq. (1^ . 


APPENDIX B: SOURCE AT A FINITE 
DISTANCE 

Here we consider a point source at a finite distance. 
For a point source, 

, S(r - To) 5{9 - n)5(z)e-^-\ (Bl) 
(I - AGfi)ro 

where (I — 4G/r)ro = y/—g. We consider a solution writ¬ 
ten in the form of the following expansion, 

OO pQQ 

(j)(r,9,z) = '^ / dk fjn,k(r) cosm9 . (B2) 

The solution for fk,m{T) is the same as fmir) in 0 but 

oj contained in ^ and are here replaced with 

and 


Nrr 


1 e^(-l)^ 

I — A Sin 


(B3) 
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OO 

First we compute (j>k{T, 9) /m,fc(?') cosm6 for r < 

m 

To- As in the case of Bessel function, we also use the 
integral representation for Hankel function 




1 

ZTT 


dse 


z sinh s—vs 


/ Ch 


(B4) 


1/Co- Performing gaussian integral, we obtain 


/'fc(r, 9) 


167r3(l - A) ^ ^-f(t',so) 

^ (l g-T^(«'+*^/2)+*e“ + ^ j • 

(BIO) 


Here the integration is to be performed along the path 
Ch presented in Fig. ^ Using the above formula and 
(112), we have 


°° ■ /•oo+iTV 

cj,,(r,9) = Y, 


327r3(l - A) 

^00+27r 


/ OO+ITT 

-OO 


' —OO 
2 


^00+27r 


167r3(l - A) 

^OO+ZTT < 

X / 






+ (0^-0) . (B5) 


We introduce a new variable t' = t + s — i7r/2. Under the 
assumption that Co ^ 1, the integration over s is domi¬ 
nated by the contribution around s = in 12. Hence, the 
integration contour for t is unaltered even if we change 
the integration variable from t iot' . After this change of 
the variable, we have 


(l>k{r,9) 


^OO+ZTT 


nOO + i-K 


167r3(l — A) 

X 


dt' 


' OO — ITT 

1 


1 _|_ g-TdA(‘'+“/2)+i®-' 


+ {9^ -9) 


(B6) 


where 


f{t', s) := C sinh(t' — s -|- i7r/2) + Co sinhs. (B7) 


We expand the exponent around a zero of its derivative. 
The derivative vanishes at s = sq, and sq is given by 


tanhso = 


C cosh(t' -|- z7r/2) — Co 
C sinh(t' -I- i7r/2) 


(B8) 


Taylor expansion of f{t', s) around s = sq becomes 


Further, we expand f{t', Sq) around an approximate sta¬ 
tionary point at t' = in 12. Then we have 

f{t',so) = *{C + Co)+ 2(^^^ ) ~ *' 2 ) (Bll) 

Again we truncate this expansion at the quadratic or¬ 
der for the same reason as before. Then one finds that 
ij)k {r, 9) is approximately given by 


==</fc(r,0), 

Stt^VC + Co 


(B12) 


where 


and 


^fc(r, 9) = , (B13) 


T{x,9) := — 


1 


2*7r(l-A) 

1 

X 


OO+ZTT 

di! 


1 — e 1 - 


(t'+i7r/2)+i& — e 


+ (9^- 


V' 

(B14) 


The function (j){r, 9) is almost identical to (j){r, 9) dis¬ 
cussed in Sec.HI, except that and other C are replaced 
with and respectively. 

Finally, we perform the integration over k. From 
(EH, we have 


C + ^o / 

(B15) 

Since C + Co = — k'^{r + ro), we can invoke the saddle 

point method again to perform fc-integral when r -|- is 
large. Evaluating the contribution from the saddle point 
at fc = Loz/D with D = ^+ {r + we obtain 


(j){r,9,z) « - 


Stt^ 


dk 


gi(«-H«o) 

VC + Co 


Akz 


T 


/(t',s) = i^ C^ + Co - 2CCo cosh(t' -k in 12) 

X ^1-I- - (s — So) -|-•••^. (B9) 

We truncate this expansion at the quadratic order be¬ 
cause the higher order terms are suppressed by 1/C or 


The calculation for r > Tq can be done in a completely 
parallel way, and the final result becomes identical to the 
case with r < Tq. 
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